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A useful choice of gauge when including null infinity in the computational domain is scri-fixing, 
that is, fixing the spatial coordinate location of null infinity. This choice allows us to avoid the in- 
troduction of artificial timelike outer boundaries in numerical calculations. We construct manifestly 
stationary scri-fixing coordinates explicitly on Minkowski, Schwarzschild and Kerr spacetimes. 

PACS numbers: 04.20.Dm, 04.25.Ha, 04.20.Jb 



I. MOTIVATION 

In numerical studies of the initial value problem for test fields on asymptotically flat background spacetimes, one 
typically truncates the solution domain by introducing an artificial timelike outer boundary into the spacetime. The 
solution then is calculated on a finite, spatially compact domain. The boundary of this domain is not part of the 
physical problem. Therefore, one tries to construct boundary conditions that correspond to transparency of this 
artificial outer boundary. In addition, these boundary conditions are required to form a well posed initial boundary 
value problem (IBVP). 

In general, it is not possible to construct such boundary conditions. Spurious reflections occur from the outer 
boundary even for the simple case of the flat wave equation on a three-dimensional ball [5(|. One therefore tries 
to minimize the amount of such spurious reflections in a manner that also ensures the wellposedness of the IBVP. 
Such boundary conditions are called non-reflecting or absorbing. It is, however, very difficult to construct them when 
the curvature of the background does not vanish or when non-linear terms appear in the equations. One needs to 
account for backscatter off curvature or self-interaction of the field near the boundary. A bad choice of boundary 
data can destroy relevant features of the solution. It has been shown, for example, that a certain choice of boundary 
data, commonly used in numerical relativity, destroys the polynomial tail behavior of solutions to wave equations on 
a Schwarzschild spacetime [8] . 

A further difficulty is related to the notion of radiation. In the study of radiative phenomena, one is interested in 
the behavior of solutions in the far-field zone. The electromagnetic field on flat space, for example, admits a global 
separation into near field terms and far field terms which simplifies the discussion of radiation considerably. When 
the background curvature does not vanish, such a separation can only be expected in the asymptotic limit to infinity 
in null directions. Therefore to study radiation accurately, one needs to calculate very large portions of spacetime. 
Using Cauchy-type foliations and truncating the computational domain, however, does not allow us to calculate such 
large portions of spacetime that astrophysically realistic distances can be modeled in a numerical calculation. 

Additional complications arise when one tries to construct absorbing boundary conditions for the Einstein equations 
@> OS BUj HH • The treatment of the asymptotic region in numerical relativity is an important theoretical and practical 
problem as it effects the accuracy of the calculated radiation field [U, |4t| . With increasing sensitivity of detectors 
and computer power, systematic errors in current numerical simulations due to wave extraction may restrict the 
predictable power of numerical relativity and blur the comparison between observed and calculated waveforms 0, |30| . 

A clean solution to the above mentioned difficulties is available on the level of geometry. The solution is to include 
null infinity in the computational domain. This geometric idea goes back to Penrose and has been known for a long 
time [42l |43|. Its implementation within the 3+1 approach, however, has been exceptionally difficult, numerically as 
well as analytically. In this article we discuss in detail the construction of suitable explicit gauges on asymptotically 
flat background spacetimes such that the geometric idea of null infinity can be incorporated in numerical calculations 
within the 3+1 approach. While we assume the spacetime to be given, a basic motivation is to extend the suggested 
methods to the treatment of the asymptotic region in numerical solutions to the Einstein equations. 

The article is organized as follows: In section [II] we introduce the basic concepts and explain the idea of scri-fixing. 
In section IIIII we present an explicit construction of scri-fixing coordinates in spherical symmetry. This method is 
applied to Minkowski spacetime in section [TV] and to extended Schwarzschild spacetime in section [V] Going beyond 
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spherical symmetry, we present the construction of a scri-fixing gauge on Kerr spacetime in section IVII We conclude 
in section IYIII with a discussion of scri-fixing and an outlook for future work. 

II. THE GENERAL IDEA 

If a background spacetime (M, g) is asymptotically flat in a suitable sense, one can attach to the spacetime manifold 
M. a boundary in null directions, such that a certain conformal rescaling of the physical metric g results in a 
metric g = fl 2 g with a smooth extension through that boundary and the conformal factor fl has a specified behavior 
on that boundary. In particular, we require > 0, Q\y = and d£L\j; ^ 0. The extended manifold with boundary 

is denoted by M.. We call the triple (A4,g, fi) a conformal extension and say that the spacetime (A4,g) is weakly 
asymptotically simple. All known explicit solutions to the Einstein equations that may be considered as "physically 
reasonable" and asymptotically flat are weakly asymptotically simple [3l| . 

By construction, null rays reach J? for an infinite value of the affine parameter along them. Therefore we refer to 
J* as null infinity. It turns out, under reasonable conditions, that J? consists of two parts ^ + and Jf~ , called future 
and past null infinity respectively, each with topology ix S 2 [43| . We will mainly focus on future null infinity. 

The propagation of massless fields is governed by the null cone structure of spacetime. It is therefore convenient in 
studies of radiative properties of test fields to choose a foliation that is closely related to outgoing null directions. An 
option is to choose a characteristic foliation. This approach has been very successful in situations where null foliations 
cover the computational domain in a regular fashion 0, I2H , |25l l36l . [53| . The main difficulty with this approach can be 
traced back to the fact that characteristic foliations are not well-behaved in regions of strong dynamical gravitational 
fields due to formation of caustics in the bundles of light rays generating the null hypersurfaces [2(| • 

Spacelike foliations are more flexible than characteristic ones. It turns out that spacelike surfaces can be constructed 
that extend through null infinity. Such surfaces are called hyperboloidal as their asymptotic behavior is similar to 
standard hyperboloids in Minkowski spacetime Instead of approaching spatial infinity as Cauchy surfaces do, they 
reach null infinity, J^, which makes them suitable for radiation extraction. It should be emphasized that hyperboloidal 
surfaces extend through null infinity as spacelike surfaces. It would be misleading to refer to them as asymptotically 
null because their tangent vector never becomes null. The term asymptotically null may be reserved for surfaces that 
become null at or close to null infinity, in which case, for example, their mean extrinsic curvature would be unbounded. 

One can formulate a Cauchy problem based on hyperboloidal surfaces [13, HI] • Although the hyperboloidal initial 
value problem and its merits have been known for a long time, numerical studies of test fields based on hyperboloidal 
foliations are rare. One major difficulty is related to the choice of gauge. Future null infinity is an ingoing null surface. 
This property can be made manifest in the gauge such that null generators of J? + converge. Such coordinates can be 
useful in numerical studies for discussing global properties of spacetimes, for example with respect to the investigation 
of a point corresponding to timelike infinity i + [33, 55]. These coordinates, however, do not seem to be convenient for 
calculating radiation along as they imply a resolution loss in the physical part of the conformal extension and 
require a gauge and boundarytreatmcnt in the unphysical part. We emphasize that, as also follows from abstract, 
semi-global existence results [la ], the resolution loss is not an intrinsic deficiency of conformal compactification, but 
a consequence of a bad choice of gauge for the type of problem that one is interested in. 

A suitable class of gauges on hyperboloidal foliations for studying radiation is given by scri-fixing 16]. In such a 
gauge the spatial coordinate location of null infinity is independent of the time coordinate so that is represented 
by the same grid points during time evolution and no resolution loss appears. This allows us to calculate radiation 
emitted to future null infinity in a very accurate way using a simple extraction procedure along an a priori known 
grid surface. Further, one can let the numerical outer boundary coincide with .y + so that no boundary treatment 
on the continuum level is required [35| . It should be noted that in a characteristic approach based on a double null 
gauge, scri-fixing is automatically satisfied by the coordinates. This seems to be one of the key advantages in including 
^ + in the computational domain within the characteristic approach. Within the hyperboloidal approach, the larger 
freedom of gauge choices needs to be carefully dealt with to achieve scri-fixing in the asymptotic domain. 

The first explicit construction of a scri-fixing gauge has been given for Minkowski spacetime [34 , 39f . Its usefulness 
has been demonstrated in studies of dynamical magnetic monopoles on a flat background [lj, Il5| . The explicit 
constructions presented in this article should be useful for studies of test fields on black hole spacetimes [54[. They 
are also intended to provide examples for testing ideas on treating the hyperboloidal initial value problem for the 
Einstein equations as well as for constructing hyperboloidal initial data [J [5a ] . 

A convenient way to discuss global properties of foliations is to depict their embedding in Penrose diagrams. We 
will employ such diagrams in this article to plot results of calculations. In spherical symmetry, Penrose diagrams 
accurately represent the global causal structure of the spacetime under study (see Q for a rigorous discussion of 
Penrose diagrams). 
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III. SCRI-FIXING IN SPHERICAL SYMMETRY 

The main property of a scri-fixing gauge is that the zero set of the conformal factor in terms of local coordinates is 
independent of time. There are two steps involved in the construction of such a gauge: the choice of a hyperboloidal 
foliation and a conformal compactification, that is, the introduction of compactifying spatial coordinates in combina- 
tion with a suitable conformal factor. We present in the following how these steps can be carried out explicitly for 
asymptotically flat, spherically symmetric spacetimes. The techniques can be easily carried over to spacetimes with 
less symmetries as demonstrated in section [VTl 

A. Hyperboloidal foliations 

Assume that we are given a globally hyperbolic, asymptotically flat, spherically symmetric spacetime (M,g) so that 
the group 5*0(3) acts non-trivially by isometry on (M,g). We can introduce on the quotient space Q = A4/SO(3) a 
natural function r : Q — > R, called the area-radius, such that the group orbits of points p £ Q have area 47rf(p) 2 . The 
metric g on A4 can then be written as g = q + f 2 da 2 , where da 2 is the standard metric on S 2 and q is a Lorentzian 
metric of rank 2 that descends naturally to a metric on Q. Due to the Jebsen-Birkhoff theorem (3?j], we can assume 
that, beside the 5*0(3) generators, there is an additional hypersurface orthogonal Killing vector held on M. that is 
timclike in the exterior [ll| . We can introduce a time function t such that this Killing vector field is given by d; and 
level sets of i are orthogonal to it. In these coordinates, the spacetime metric reads 

g = qn dP + qff df 2 + f 2 da 2 , ( 1 ) 

where qn and qff are functions of f only. The time function t is by construction such that the time slices are Cauchy 
surfaces and approach spatial infinity in the asymptotic region. We are interested in null infinity. Therefore, we need 
to introduce a new coordinate t = t(t,r) such that the time slices are hyperboloidal surfaces. 

In stationary spacetimes one may define the notion of Killing observers. The worldline of a Killing observer 
corresponds to the integral curves of the timelike Killing field. These observers are not inertial if the curvature of the 
spacetime does not vanish, but they are natural in the sense that their distance to each other remains constant. A 
Killing observer far away from a source can be regarded as representing a gravitational wave detector. Therefore, we 
wish to keep the natural coordinate representation of the Killing observers so that the Killing vector field is given by 
dt in terms of the new coordinate t. This implies a transformation of the form 

t = t-h(r), (2) 

where h(f) is called the height function. In this new time coordinate, the quotient metric q can be written as 

q = (Judt 2 + 2 j b! dt df + (q^ h' 2 + qff. ) df 2 , where b! := — . 

All metric components are functions of f only. A major advantage of this manifest time-independence of the metric 
functions is the implied simplification in numerical applications. 

The next question is how to choose the height function to obtain a hyperboloidal foliation. As the time slices are 
required to be spacelike and to approach future null infinity, the height function should lead to a " bending up" of the 
surfaces as f — > oo such that, asymptotically, their behavior is similar to outgoing null surfaces. In the next subsection 
we will derive a necessary requirement on the asymptotic behavior of the height function for the regularity of the 
conformal compactification. 

B. Scri-fixing conformal compactification 

On level sets of the hyperboloidal time function t we introduce a compactifying coordinate r. The physical metric 
g then becomes singular at the domain boundary of the compactifying coordinate. To compensate for this singular 
behavior we rescale the metric with a conformal factor Q that vanishes at that boundary in a suitable way. This 
procedure is referred to as conformal compactification. 

It is natural and convenient to require that the compactifying coordinate r is an area radius on the conformal 
extension with respect to the rescaled metric g = fl 2 g. This requirement implies f = r/fi. The zero set of the 
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conformal factor corresponds to the asymptotic region of the physical spacetime. Writing the rescaled quotient metric 
q = n 2 q in terms of the compactifying coordinate r using r(r) = r/i7(r), we get 



q : ■ V/r lu <// 2 - 2 (/,-, //' (<> - <>' /■)<//,/;- -I- y~tt h ^~y ff (J) . {)' r r,jr 



This formula can be applied with almost any choice of conformal factor. Note that we have some freedom here. We 
can require, for example, that r be the proper distance in the conformal extension, i.e. q rr = 1. Then the conformal 
factor is determined by this requirement and cannot, in general, be written in explicit form [40| . In contrast, by 
allowing the metric component q rr to be more general, we are free to prescribe the representation of a conformal 
factor in terms of a suitable compactifying radial coordinate r and the coordinate transformation is explicit. 

A convenient representation for the conformal factor that we will resort to is O = (1 — r). Our discussion is not 
restricted by this choice. Other coordinate representations of the conformal factor are also allowed. The above choice 
is not good at the origin as the gradient of r vanishes at this point. A suitable representation can be used in the 
interior as will be discussed in the next subsection. The above choice fixes the compactifying coordinate r via 



implying 



r(r) = — (3) 



q = n 2 q s dt 2 + 2 q E ti dt dr + dr 2 . (4) 



The vanishing of qtt at {fi = 0} implies that dt becomes null at this surface. For regularity of the metric we require 



h' = t — — as f — > oo, with C = const. (5) 

V lit 

To gain an intuition for the behavior of hyperboloidal surfaces, we will calculate coordinate speeds of characteristics. 
The coordinate speed c of ingoing and outgoing null rays in radial directions can be calculated by c T = {—qtr ± 
\J Qtr ~ Qrr qtt) / qrr- The speeds are independent of the conformal rescaling of the metric which is a reflection of 
the fact that conformal rescaling leaves the null cone structure invariant. The positive sign on the right hand side 
corresponds to ingoing null rays, the negative sign corresponds to outgoing null rays. On the conformal boundary we 
get c T = (—qtr ± qtr)/qrr- As expected, the ingoing characteristic speed is zero. Consequently no analytic boundary 
treatment is needed or allowed. The outgoing characteristic speed is c+\y+ — —2q tr /q r r- In a future hyperboloidal 
foliation we have qtr\,<f+ < and q rr > 0, so that c+|j?+ > and radiation can leave the spacetime through future 
null infinity. 



C. Matching 

An important feature of hyperboloidal surfaces, in contrast to characteristic ones, is that their distinctive feature is 
determined in the asymptotic region. While being spacelike or null is a local property, being hyperboloidal or Cauchy 
is a global property. This fact allows us to construct useful hyperboloidal foliations easily by a matching procedure. 

Foliations commonly used in the 3+1 approach are characterized by properties which imply that the hypersurfaces 
approach spatial infinity in the asymptotic domain. One may wish to keep these properties in the interior and 
transform the surfaces in the exterior such that they approach null infinity asymptotically. Hyperboloidal foliations 
can be constructed that coincide with arbitrary spacelike foliations in the interior. To achieve this, we will employ 
the following smooth transition between interior and exterior height functions denoted by hi and h e respectively 

)hi(f) for f<fi, 

e-(?-n) 2 /(f ; -' ; e ) 2 (^(f)_/ le (f)) + / le (r) for fi<f<r e , (6) 
h e (f) for f >f a . 

The height function needs to satisfy \h'\ < \/qff /(—qii) f° r the foliation to be spacelike. We may apply the above 
transition also with respect to derivatives of the corresponding height functions. We may set the conformal factor to 
unity in the interior with a smooth transition to its asymptotic representation in the exterior. The transition domain 
(fj,f e ) should be in the exterior, but not necessarily in the far-field zone. The free parameters of the matching can 
be decided upon by numerical experiments. Detailed examples will be discussed in later sections. 

One should note that this matching is fundamentally different from Cauchy-characteristic matching where one 
matches a spacelike and a null surface. The matching in our case results in a single, smooth, hyperboloidal surface 
with an arbitrary behavior in the interior determined by hi(f). 
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IV. MINKOWSKI SPACETIME 



A. Cauchy foliation 



The 2-dimensional Minkowski metric fj on the quotient space in standard coordinates is given by f/ — —dt 2 + df 2 . 
The Penrose diagram in Fig. [1] depicts the causal structure of Minkowski spacetime. Each point in the diagram 
represents a sphere except the point i that corresponds to spatial infinity and the points along the vertical line 
segment connecting i + and i~ that correspond to the origin {f = 0}. Radial null surfaces are represented by straight 
line segments with 45 degrees to the horizontal. Curves with a smaller angle represent spacelike surfaces, curves with 
a larger angle represent timelike surfaces. Level surfaces of f approach timelike infinity denoted by V s1 , level surfaces 
of t approach spatial infinity denoted by i . 





Fig. 1: Penrose diagram of Minkowski spacetime. Solid 
lines represent t — const, surfaces, dashed lines repre- 
sent f — const, surfaces 



Fig. 2: Characteristics in Minkowski spacetime in the standard foli- 
ation with respect to the compactifying coordinate r defined in ([3| . 
In- and outgoing null surfaces are represented by dashed and solid 
lines respectively. Note the different scales in t and r. 



It is evident from Fig. Q]that an outgoing null ray never leaves the surfaces of a Cauchy foliation. This is demon- 
strated more clearly by the characteristics with respect to coordinates of the standard foliation. We introduce on the 
Cauchy surfaces the compactifying coordinate r defined in ([3]). The resulting metric reads fj = —dt 2 + dr 2 /(l — r) 4 , 
and the coordinate speeds of characteristics become c± = ±(1 — r) 2 . The characteristic speeds decrease to zero in the 
asymptotic region close to r = 1. In simulations, this results in a slowing down of outgoing waves. This feature is 
clearly shown in Fig. [2] where in- and outgoing characteristics are depicted by dashed and solid lines respectively. We 
see that outgoing characteristics pile up during the evolution near the boundary and eventually, cannot be resolved 
(see [45[ for a discussion of this phenomena in the context of a binary black hole evolution) . 

A suitable compactified foliation should not only have no ingoing characteristics coming in from the outer boundary, 
but also allow outgoing characteristics to leave the simulation domain. 

B. Hyperboloidal foliations 

A typical example for a future hyperboloidal foliation is plotted in Penrose diagram Fig. [3J Null rays do leave the 
surfaces of a hyperboloidal foliation. The corresponding time function to the foliation in Fig. [3] is given for t > by 



f 2 



i 



(7) 



A distinctive feature of hyperboloidal surfaces is that their mean extrinsic curvature attains a non-vanishing, finite 
value at infinity. We use the convention in which positive K corresponds to expansion so that positive asymptotic 
values of K correspond to future hyperboloidal surfaces. The mean extrinsic curvature of level sets of t as defined 
in ([7]) reads K — 3/vI + 1 2 . Each surface has a spatially constant mean curvature. The foliation, however, does 
not seem to be convenient for numerical calculations, because the time function does not have the form ([2]). The 
Minkowski metric in this new coordinate depends on time. 

A more convenient foliation can be constructed by Lie dragging a standard hyperboloid along the timelike Killing 
field. The spherically symmetric standard hyperboloids in Minkowski spacetime are given by a 2 = t 2 — f 2 , with 
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Fig. 3: Penrose diagram of a hyperboloidal foliation of Minkowski spacetime. Solid lines represent hyperboloidal surfaces, 
dashed lines represent Killing observers 



a £ R+. Translating these surfaces along by an amount of t, we get the equation a 2 = (t — t) 2 — f 2 . The new time 
coordinate t is then defined by 



t(t,f) = i- \Ja 2 + r 2 



(8) 



Each surface of the foliation has a constant mean curvature K = 3/ a. The Penrose diagram in Fig. [4] shows foliations 
defined by {5J for different values of K. 






Fig. 4: Penrose diagrams for constant mean curvature (CMC) foliations of Minkowski spacetime. Solid lines represent hyper- 
boloidal surfaces, dashed lines represent Killing observers. The foliations have K = {6, 3, 2} from left to right. 



In this new time coordinate, scri-fixing conformal compactification as in results in 



V 



-n 2 dt 2 



2r 



■. dtdr 



a 2 n 2 



dr 2 



(9) 



\]a 2 tt 2 + r 2 

We see from the Penrose diagram Fig. [4] that the mean extrinsic curvature of the surfaces determines the angle at 
which the surfaces cut null infinity. It also acts as a measure for how close the surfaces are to outgoing null surfaces. 
The similarity of future hyperboloidal surfaces to outgoing null surfaces can be quantified by the coordinate speed 
of the outgoing characteristic at null infinity. We have c+|j?+ = 4K 2 /9. A large value of K corresponds to surfaces 
which are closer to null surfaces, as also suggested by the Penrose diagrams in Fig. [4] This is seen more clearly in 
Fig. [5j where in- and outgoing characteristics have been plotted with respect to coordinates. Note that the choice of a 
high value for if in a numerical calculation using an explicit time integration algorithm would require very small time 
steps due to the Courant condition which demands that the numerical cone needs to be larger than the characteristic 
cone. For small values of K, however, this condition can be easily satisfied as seen also in Fig.O 

It should be emphasized that the notion of " proximity" to null surfaces does not have a general geometric meaning. 
Geometrically, it makes no sense to claim that a spacelike surface is close to a null surface as this depends on the 
chosen frame of reference. On the other hand, in stationary spacetimes, we have a natural choice of observers, namely 
the Killing observers, so that we can introduce such a notion. In the general non-stationary case, this notion may be 
useful in the asymptotic domain with respect to the asymptotic Killing observers. 
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Fig. 5: Coordinate representations of in- and outgoing characteristics of with K = 3/a = {6, 3, 0.5}. High values for K are 
seen to be restrictive in a numerical calculation in which the numerical cone needs to be larger than the characteristic cone. 



C. Matching 



Level surfaces of the standard time coordinate t in Minkowski spacetime are maximal, i.e. their mean extrinsic 
curvature vanishes. To match an interior maximal foliation to a hyperboloidal foliation, we choose in j6]) 

hi = d, h e = C e + ^/a 2 + P. (10) 

The constants Cj and C e must be chosen such that the matched surface is spacelike in the transition domain, that 
is, \h'\ < 1. One convenient way to achieve this is to require that the surfaces intersect within the transition domain. 
We denote this matching point by r m € (fj,f e ), and require hi(f m ) — h e {f m ). This implies C e = C% — y/a 2 + f^,. 
The resulting foliation has been plotted in Fig. [6] The parameters have been chosen such that the matching is clearly 
visible in the diagrams. In a numerical calculation, they need to be decided upon by numerical experience. 




Fig. 6: Transition between maximal surfaces and two types of hyperboloidal surfaces determined by (|10|l and {jXTJ . The three 
dashed lines depict the timelike surfaces fj = 1/2, r m = 1, and r e — 2. The exterior foliation has a = 2 in both cases. 

A more general way of constructing hyperboloidal surfaces in a given spacetime is to Taylor expand the formula 
([5]) near infinity. We get in Minkowski spacetime 

/ a 2 \ 3 a 2 a 2 / 1 \ 

^-(^-(TTFj*) =1 -^P + P +0 {^)> as ^°°> ( n ) 

where we set C = a 2 . The first two terms of the expansion are sufficient to fulfill the requirements on the height 
function for a future hyperboloidal foliation. The foliation is defined in the exterior domain for f > a — 1, or written 
in the compactifying coordinate, for r > 1 — 1/a. We use ((6]) to match the surfaces to an interior maximal foliation 
by choosing h e (f) = r + a 2 / (2f) + C e , and h{ = The resulting foliation has been plotted in Fig. [S] The rescaled 
metric in the exterior domain reads 

1] = -n 2 dt 2 -2(1-^-) dtdr + irdr 2 . 

The outgoing characteristic speed at future null infinity is c + \y+ — A /a 2 . The mean extrinsic curvature of these 
surfaces depends on r and attains the value K = 3/a at infinity. 
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V. SCHWARZSCHILD SPACETIME 



A. Cauchy-type foliations 

All foliations of Schwarzschild spacetime commonly employed in numerical calculations within the 3+1 approach 
are Cauchy-type foliations approaching spatial infinity in the asymptotic domain. We mention some of them in this 
subsection for later reference. We focus on the physically interesting part of the Schwarzschild-Kruskal extension (the 
regions I and II in [311 ]) which we refer to as the extended Schwarzschild spacetime. 



1, Schwarzschild coordinates 



The 2-dimensional Schwarzschild metric in Schwarzschild coordinates is given by 



2m\ ,~ 2 / 2m \ _> 



q = -\l-^-jdf+\l-—\ df, f>2m. (12) 

Level surfaces of t are maximal. They intersect at the bifurcation point and at spatial infinity as seen in Fig. [7] 
Therefore, these coordinates are suitable neither to study the region near the black hole, nor the asymptotic domain. 



2. Eddington-Finkelstein coordinates 



In numerical applications, a common choice of so-called "horizon penetrating" coordinates in Schwarzschild space- 
time is defined as in ^ with the height function 



h(f) = 2m In 



r 

2m 



The coordinates such defined have originally been discussed by Eddington [T(| and Finkelstein [l3| . The Schwarzschild 
metric on the extended Schwarzschild manifold takes the form 



2m 
r 



dt 2 



4m 



■ dtdf 



2m 

f 



dr 



,-.2 



f > 0. 



(13) 



The mean extrinsic curvature of level sets of t depends on f and vanishes in the asymptotic domain. The foliation by 
the time function t is depicted in the Penrose diagram Fig. [H We see that the foliation penetrates the horizon which 
makes it suitable for excising the singularity from the computational domain. The surfaces are spacelike at each point. 
They approach spatial infinity in the asymptotic domain which makes them inconvenient to study radiation. 



singularity 




Fig. 7: Level sets of Schwarzschild coor- 
dinates on the extended Schwarzschild 
manifold. 



Fig. 8: Eddington-Finkelstein folia- 
tion suitable for the excision tech- 
nique. Dashed lines correspond to 
r = {1.5m, 4m}. 



Fig. 9: Maximal slices suitable for the 
moving puncture technique. Dashed 
lines correspond to f — {1.5m, 4m}. 
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3. Maximal slices for the puncture technique 

Currently, the most popular approach in numerical calculations of black hole spacetimes is the moving puncture 
technique. The evolution with puncture initial data for a Schwarzschild black hole in a modified 1+log slicing reaches 
a final stationary state that corresponds to a maximal slice given in suitable coordinates 0, H3, HH 0] • We will use 
maximal foliations to demonstrate the possibility that hyperboloidal surfaces can be constructed which are suitable 
for the application of the moving puncture technique in the interior. 

The requirement that the mean extrinsic curvature K of the surfaces defined as in ([2]) are constant implies in 
Schwarzschild spacctime 




(14) 



(l — 222.) 1 — h' 2 (l — — ) 

We can integrate this relation once for a maximal foliation (K = 0) in Schwarzschild spacetime to obtain 0, [HJ 

C 

h ' {f) = "(1- 2 f )V^--2mf3 + C 2 ' (15) 
where C is a constant of integration. The Schwarzschild metric becomes 



/ 2m\ , 2 2C , ,„ 1 ,_, . , , / 2m C 2 

q = -[l-— J dt 2 + — dtdr + -dr 2 , with f = ] jl- — + —. 

The standard foliation of Schwarzschild spacetime is a member of this family with C = 0. The final stationary state 
for a puncture evolution in Schwarzschild spacetime is similar to a limiting maximal slice with C = 3V3m 2 /4. Fig. [9] 
shows a maximal foliation with C — 3\/3m 2 /4 in the Penrose diagram of the extended Schwarzschild spacetime 
together with the limiting surface at f = 3m/2 where / = 0. 



B. Hyperboloidal CMC-foliations 

The main difference between the Schwarzschild foliation in Fig. [7] and the horizon penetrating foliations in Fig. [8] 
and Fig.[9]is the behavior of surfaces near the event horizon. The intersection at the bifurcation point in Schwarzschild 
coordinates is not suitable to study the black hole. Similarly, the intersection of surfaces of a foliation at spatial infinity 
is not suitable to study radiation near null infinity. The construction of hyperboloidal foliations, while geometrically 
quite different, can in this respect be regarded as similar to the construction of horizon penetrating coordinates. 

To find a representation of Schwarzschild spacetime based on hyperboloidal surfaces, we use the general family of 
spherically symmetric CMC surfaces in the extended Schwarzschild spacetime as discussed in [(| [38| (see [5l| for a 
numerical study). The integration of (fTl)) with K ^ and integration constant C leads to the differential equation 



ti 



K r_ 

3 



c 



(1 



P(f) 



with P(f) 



( kf 3 



c 



2m 



(16) 



The maximal foliation given in (jT5j) is a member of this family with K = 0. To get the height function we need to 
integrate the differential equation in Q16p . The integral cannot be written in explicit form but needs to be carried out 
numerically. The resulting foliation has been plotted in Fig. [TDJ 

Scri-fixing conformal compactification as in ^ based on this foliation results in 



1 = ~ 1- 



> n 2 dt 2 - v 



P(r) 



dtdr 



P 2 (r 



■dr 2 , 



(17) 



where 



P(r) :=Q 3 P(r) = A (*lt-CW 



2m(l 



Q 2 r 4 
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Fig. 10: CMC-foliation of the extended Schwarzschild spacetime with K — 0.4, C — 2 and m — 1. On the left, the Penrose 
diagram is plotted, on the right in- and outgoing characteristics. The excision surface is taken to be at f = 1.8m and is plotted 
as the dashed line in the Penrose diagram. The event horizon is at r — 2/3. 

The global behavior of CMC-surfaces in Schwarzschild spacetime depends on the free parameters as discussed in 
[38| . The Penrose diagram in Fig. [10] depicts surfaces which are convenient for numerical calculations. They come 
from future null infinity, pass the horizon above the bifurcation sphere and run into the future singularity. Depending 
on the free parameters, P(r) may have a non-vanishing zero set. The corresponding surface lies always inside the 
event horizon and does not pose a difficulty if the excision surface has been chosen carefully. 

The outgoing characteristic speed at ^ + reads c+\j?+ = AK 2 /9 : just as for the CMC-foliation of Minkowski 
spacetime. The main difference to Minkowski spacetime is the black hole region. The coordinate representation of 
characteristics in Fig. [10] illustrates that there are no incoming characteristics into the simulation domain both at the 
excision boundary inside the event horizon and at the outer boundary on future null infinity. 

C. Matching 

Schwarzschild spacetime is asymptotically flat in null directions, implying that its asymptotic structure is similar to 
Minkowski spacetime. We can therefore use the same form of the height function as we did for constructing a hyper- 
boloidal foliation in Minkowski spacetime to get, at least asymptotically, a hyperboloidal foliation in Schwarzschild 
spacetime. In writing the height function, however, we must be careful about the choice of coordinates. Asymptotic 
flatness in null directions can be made manifest in terms of the Reggc- Wheeler coordinate f„ = r + 2m ln(r — 2m) 
as described in (44|. It is therefore the height function h(r) = a 2 + fj, that gives a hyperboloidal foliation in the 
exterior domain. As the coordinate f* is not suitable near the black hole region, we match the surfaces to interior 
Eddington-Finkelstein surfaces. The resulting foliation has been plotted in Fig. [TT1 




Fig. 11: Penrose diagrams of matched foliations in the extended Schwarzschild spacetime. The dashed lines correspond to 
f = 1.5 m, fi = 2.1m, f m = 2.5 m, r e — 3 m. The diagram on the left illustrates matching between interior Eddington- 
Finkelstein surfaces and surfaces with height function h = The diagram on the right illustrates matching between 
maximal slices and hyperboloidal surfaces with the height function (|19l) . 
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The simplest method to construct hyperboloidal surfaces seems to be the Taylor expansion of (0 as in (fTTj) . In 
Schwarzschild spacetime, the expansion delivers 

L ,,_. f / / 2m\ C \* l 2m 8m 2 - C C + 8m 3 

h(r)~- — —l-l-_- — . =i + _ + — — + + o asr^oo. 18) 

w r-2m \ V r J (1 + r) 2 / r 2r 2 r 3 \ r / 

We take the first two terms so that the height function becomes 

h(r) = f + 2m\nf, (19) 

which essentially corresponds to the Regge- Wheeler coordinate in the exterior domain. Matching to an interior 
maximal foliation gives the second conformal diagram in Fig. 1111 The conformally rescaled metric in a scri-fixing 
gauge in the exterior domain reads 



2mfl\ „ 2 , o / Am 2 Q 2 \ , , 8m 2 (r 2 - 2m 2 Q 2 ) , „ 
-—) V-dt 2 -2(l-— ytdr + r3 \ r _ 2mQ) > dr 2 . 

The outgoing characteristic speed at future null infinity is c+\ji+ = l/(2m 2 ). The mean extrinsic curvature of the 
time slices attains a positive finite value at infinity that depends on the mass of the Schwarzschild spacetime. It reads 
K\jr+ = 3/(2V2m). To gain some control on properties of this asymptotic foliation one can take further terms in the 
Taylor expansion (fTS)) . 



VI. KERR SPACETIME 



In most studies, going from the spherically symmetric case to the axially symmetric case is a challenging task. Take 
as an example the construction of suitable null coordinates. Null coordinates can easily be constructed on Minkowski 
and Schwarzschild s pac etimes. In Kerr spacetime, however, they are difficult to deal with and cannot be written 
in explicit form [4rl |52|. This is, in part, a consequence of the fact that null coordinates are very rigid. Spacelike 
foliations on the other hand are simpler to construct due to the large freedom involved in their choices. Further, 
the property of a spacelike surface to be hyperboloidal does not restrict its form in the interior as we have seen in 
previous sections. It only determines its asymptotic behavior. This fact in combination with the asymptotic flatness 
of Kerr spacetime allows us to carry over the construction of a scri-fixing gauge in spherical symmetry directly to 
axial symmetry on the example of Kerr spacetime. 

The most commonly used coordinates for Kerr spacetime in numerical relativity are the Boyer-Lindquist and the 
Kerr-Schild coordinates. In the following we present the asymptotic form of simple hyperboloidal foliations in these 
coordinates. Other choices are possible and straightforward using the techniques presented so far. The Kerr metrics 
as given in (|2ip and (|22p can be used in the exterior domain. One may employ the matching technique to obtain a 
convenient foliation also in the interior depending on the problem at hand. 



A. In Boyer-Lindquist coordinates 



The Kerr metric in Boyer-Lindquist coordinates reads 



where 



1 ^- dP =r- sin 2 i? dtdip + ^r dr 2 + S dd 2 + R 2 sin 2 d dip 2 , 

E J E A 



E = f 2 + a 2 cos 2 <&, A = f 2 + a 2 - 2m f, R = r + a + ^ . (20) 



By m we denote the ADM-mass and by a the angular momentum of Kerr spacetime. To choose a suitable height 
function for a hyperboloidal time coordinate in the exterior domain, we use the Taylor expansion of ([5]). It reads 

2m 8m 2 -C-a 2 sm 2 $ „ { l\ 

h'(r) = 1 + — + — + O _ as r -► 00. 

r 2r z \ r A I 
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We use the first two terms and set h(f) = f + 2mlnf, as for Schwarzschild spacetime in Schwarzschild coordinates. It 
is interesting that the height function for a hyperboloidal foliation does not necessarily depend on angular momentum 
but on mass. On the other hand, the form of the height function should not surprise us, as asymptotically, the Kerr 
metric in Boyer-Lindquist coordinates approaches the Schwarzschild metric in Schwarzschild coordinates. With the 
above choice of the height function, the components of the metric in the exterior domain of Kerr spacetime in a 
scri-fixing gauge take the form 

2mrfl\ ~ / 2mrft\ ( 2mfi\ 

2amrfl 3 9 „ 2amQ , ~, o „ , 

9t v = ^ siaTti, g rip = — (r + 2mO)sm^, (21) 

gm = E, g w = R 2 sin 2 

where we have set S = £7 2 E, A = f2 2 A, and R — flR. A calculation shows that the rescaled metric component g rr is 
manifestly regular at {f2 = 0} and can be written as 





8m 2 / mO\ a 2 sin 2 4m 2 2m , 9 „. 9 ^ 9 „.„ 9 ^ 9 , 

g rr = 1 + r 5- + — — (r 2 (2mr - a 2 fl) - a 2 fl cos 2 i?(2a 2 2 + r(r - 4mfi))) 

E \ r J A r 2 - rAE 

The rescaled metric becomes at future null infinity 



g \ J ?+ = -2dtdr+ {8m 2 - a 2 sin 2 i?) dr 2 + diT + siiT i? dyT . 

We see that, in these coordinates, the coordinate speed of outgoing characteristics depends on i? in a manner related 
to the angular momentum of the Kerr black hole. 



B. In Kerr-Schild coordinates 



The Kerr metric in ingoing Kerr-Schild coordinates reads 

/ 2mf\ 4mr Aamr . , „ , „ , ,„ / 2mf\ , , 
.9 = - 1 =- dr + dtdr Tzr— svcT'd dt dip - 2a snr# 1 + drdip 



E 



2mf \ 



1 + — j df 2 + E di? 2 + E 2 sin 2 

with the definitions (f2T))) . For the choice of the height function we use the Taylor expansion of ((SJ) as before and get 

, 4m 16m 2 -C -/IN 
h'(r) = 1 H — - — I —5 hO - , as r -> 00. 



Integrating the first two terms gives h(f) = r + 4m Inf. With this choice of the height function, the components of 
the resulting metric in the exterior domain of Kerr spacetime in a scri-fixing gauge take the form 

/ 2mrSl\ ^ 9 / 4mSl\ 4mf2 . „ ^ x 

gtt = - f 1 - — ^— J 9tr = - f 1 + — J + -g-(r + 2mQ), 

5«¥> = g sm 2 tf, 9r ^-a S m^ 1 + — H2 m fi) , (22) 



9rr = —5 — (r + 2mf2)(2mr — Via 2 cos 2 t9), = E, = R 2 sin 2 i9, 

r^E ' 



The rescaled metric becomes at future null infinity 

g\j?+ = -2dtdr + 16m 2 dr 2 - 2as'm 2 -ddr dtp + dd 2 + sin 2 dtp 2 . 
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VII. DISCUSSION 



It has long been suggested that using a hyperboloidal foliation in combination with conformal compactification of 
the metric to include null infinity in the computational domain would be a solution to the outer boundary problem 
as well as the radiation extraction problem within the 3+1 approach in numerical relativity (l7l . I32I ]. This method 
requires an adequate choice of coordinate and conformal gauge in the asymptotic domain. An especially suitable class 
of gauges for numerical studies on radiative properties of fields is scri-fixing [lol |35| . While it has been known how 
to construct an explicit scri-fixing gauge in Minkowski spacetime [HI HH , there was no such construction available 
for black hole spacetimes. In this article we have explicitly constructed and studied scri-fixing gauges in Minkowski, 
Schwarzschild and Kerr spacetimes. The presented constructions are in such a generality that the methods can be 
applied to any explicitly given asymptotically flat spacetime admitting a smooth conformal compactification at null 
infinity. 

We may argue that scri-fixing gauges are both natural and useful. Observers with respect to which the coordinate 
location of ,J? + is independent of time represent asymptotically stationary observers. We have seen in the examples 
that in the asymptotic limit, Killing observers become null. This is also seen by the fact that a scri-fixing gauge in a 
stationary spacetime can be constructed such that the metric is manifestly time independent. In this sense, scri-fixing 
reflects the astrophysical situation that the gravitational wave detectors are asymptotically stationary observers of 
an isolated system. It is also with respect to such observers that notions of conserved quantities can conveniently be 
defined. Choosing a scri-fixing gauge would considerably simplify the calculation of such quantities and the analysis 
of numerical solutions. Due to the manifest regularity of the rescaled metric in our construction, points at infinity 
can be treated on the same level as finite points. 

A further suitable property of the foliations we discussed is their intrinsic time asymmetry. Many calculations in 
numerical relativity are made with time symmetric data. It seems highly unlikely that such data can be regarded as 
physically reasonable. Hyperboloidal surfaces on the other hand do not allow time symmetry by their very nature 
and therefore reflect the physical situation more closely. 

That an explicit scri-fixing gauge can be useful has been demonstrated in [l4l.[l5| in Minkowski spacetime. A natural 
next step would be to extend such studies to black hole spacetimes and to other test fields. One can also study to 
what extent the outer boundary introduced in currently common numerical calculations influences the solution. In 
particular, the matching technique would allow us to make such a comparison conveniently. 

The main motivation to include null infinity in the computational domain is to have access to the correct asymptotic 
waveform. It seems, however, that the approximations commonly used in numerical relativity may be sufficiently 
accurate for many purposes. An interesting question to study in this context would be whether and to what extent 
numerical access to null infinity is relevant for astrophysically realistic predictions concerning gravitational radiation 
signals to be measured on Earth [24j . Just to give an example, it is known that the tail behavior is different at null 
infinity and at a finite distance away from the source I26IJ47I . In this case, a naive interpretation of finite distance 
extractions may even lead to a misleading error estimate |54| . Certain other features of gravitational radiation might 
also depend on the observers location in a non-trivial manner. 

One may gain some insight from the conformal compactifications presented in this article for the construction of 
hyperboloidal initial data |l| . The examples may also serve as a test ground for new ideas on treating the hyperboloidal 
initial value problem for the Einstein equations. The construction of a scri-fixing gauge in this case can be achieved 
by prescribing a suitable representation for the conformal factor in terms of coordinates in combination with an 
appropriate choice of a gauge source function in the general wave reduction of the Einstein equations [56] . Whether 
this approach can be implemented in numerical calculations remains to be seen. The explicit examples presented in 
this article may be useful for performing numerical experiments in this approach. 
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